**Volume Title** 

ASP Conference Series, Vol. **Volume Number** 
**Author** 

© **Copyright Year** Astronomical Society of the Pacific 



2D computations of g modes in fast rotating stars 

J. Ballot/'^ F. Ligmeres,''^ V. Prat,^'^ D. R. Reese^ and M. Rieutord^'^ 

^CNRS, IRAP, 14 avenue Edouard Belin, 31400 Toulouse, France 

^Universite de Toulouse, UPS-OMP, IRAP, 31400 Toulouse, France 

^Observatoire de Paris, LESIA, CNRS, Universite Pierre et Marie Curie, 
Universite Denis Diderot, 5 place J. Janssen, 92195 Meudon, France 

Abstract. We present complete 2D computations of g modes in distorted polytropic 
models of stars performed with the Two-dimensional Oscillation Program (TOP). We 
computed low-degree modes (i = 1 modes with radial order n - -1 ...-14, and 
{ - 2,3 modes with n - -1 ... -5 and -16 . . . -20) of a nonrotating model and fol- 
lowed them by slowly increasing the rotation rate up to 70 % of the Keplerian break-up 
velocity. We use these computations to determine the domain of validity of perturbative 
methods up to the 3rd order. We study the evolution of the regularities of the spectrum 
and show quantitative agreement with the traditional approximation for not too large 
values of the ratio of the rotation rate to the pulsation frequency. We also show the 
appearance of new types of modes, called "rosette" modes due to their spatial structure. 
Thanks to the ray theory for gravito-inertial waves that we developed, we can associate 
these modes with stable periodic rays. 



1. Introduction 

Rapid rotators are typic al among the main-sequence stars of intermediate or high masses 
{e.g. , Rover et al.l2007 ). Some stars exhibit oscillations due to waves trapped in the stel- 



lar interior. For some classes of pulsating stars, such as 6 Set or /3 Cep, these oscillations 
are p (pressure) modes, whereas for y Dor or Slowly Pulsating B (SPB) stars, they are g 
(gravity) modes. For the present study, we focus on g modes, i.e. low frequency modes 
driven by the buoyancy force. 

To be able to model the oscillation spectrum of rotating stars, the effects of rotation 
must be carefully considered. Rotation distorts stars under the effects of the centrifu- 
gal acceleration, and the Coriolis force can substantially modify the stellar oscillation 
properties. Complete computations are hard and delicate because it is a 2D eigenvalue 
problem and it requires 2D models of stars. To simplify the problem, the effects of 
rotation can be treated as a perturbation when the rotation rate, Q, is small relative to 
mode frequencies, cj, and to the Keplerian break-up rotation rate, Qk = (GM/RI^)^^'^ 
(M and /?eq are the mass and equatorial radius of the star, respectively). Nevertheless, 
these hypotheses are not in the least verified for rapid rotators. 

We performed complete two-dimensional computations of g modes in distorted 
polytropic models of stars. After briefly presenting our computations (Sect. |2ll, we 
show the domains of validity of perturbative approaches (Sect. |3]l, summarising results 
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presented in IB allot et al] (l201Clh . Then we discuss the evolution with rotation of the 
spectrum regularities (period spacing) and compare our results to those obtained with 
the traditional approximation (Sect.|4ll. Finally, in Sect. [51 we discuss the spatial distri- 
bution of eigenmodes, especially the appearance of rosette modes. 



2. 2D computations 

We consider fully radiative stars in this study. Since g modes are driven by the buoyancy 
force, they are evanescent in convective zones. SPB and y Dor stars have large radiative 
regions with a convective core, and a thin convective envelope for the latter. The effects 
of convective cores are not considered here. 

We approximate the equilibrium structure of rotating sta rs with self-gra vitating 
uniformly-rotating polytropes, as already done for example in iLignieres et all ([2006). 



Poly tropes are described in the co-rotating frame by the following system of equations: 

Po = Kp];"^ (1) 

"^Po = Pogo (2) 

A(Ao - 47rGp« (3) 

where p„ is the pressure, the density, ipo the gravitational potential, K the polytropic 
constant, ju the polytropic index, G the gravitational constant, and go the effective grav- 
ity, defined as 

go = -^{^o - d^s^ll), (4) 

where s the distance to the rotation axis. Since the star is distorted by the centrifugal 
force, it is not spherical and we used a suited surface-fitting spheroidal system of coor- 
dinates based on lBonazzola et all (1 19981) . This equation sy stem is numerically solved 
with the ESTER code (Evolution STEUaire en Rotation, see lRieutord et al.ll2005l) . This 



is a spectral code using Chebychev polynomials in the pseudo-radial direction, and 
spherical harmonics in the horizontal one. We computed models decomposed on spher- 
ical harmonics up to degree { = 32, which ensures a sufficient accuracy for our pur- 
poses. The polytropic index is fixed to /u = 3, thereby approximating the structure of a 
fully radiative star. 

In the co-rotating frame, the following equations govern the temporal evolution of 
small adiabatic inviscid perturbations of the equilibrium structure: 

dtp - -V-(poV) (5) 
Podt\ ^ -Vp + pgo - Pc^^ - 2po^ X V (6) 



d,p-cld,p = '"'" v- go (7) 

go 

Ai// - 47rGp (8) 

where p, p, v, and ip are the Eulerian perturbations of density, pressure, velocity, and 
gravitational potential, = Tipo/po the adiabatic sound speed, and A^^^ = go -(V In po- 
Vln Po/^i) the Brunt- Vaisala frequency. Fi denotes the first adiabatic exponent. 

We look for time-harmonic solutions of the system ©-([S]!, and we solve the ob- 
tained eig envalue problem with the two-dimensional oscillation program (TOP, see de- 



poNy^ 



tails in iReese et al.ll2006r) . The equations are projected on the spherical harmonic basis 
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Y"\ Due to the axisymmetry of the system, the projected equations are decoupled rel- 
atively to the azimuthal order m, but in contrast to the spherical nonrotating case, they 
are coupled for all degrees £ of the same parity. We compute modes with typically 20 
spherical harmonics (^max = 40 + |m|). 

We computed frequencies, a>||'^, of ^ = 1 to 3 modes in a nonrotating polytrope (in 
this case, modes are described by one degree only). We considered £ = I modes with 
n = -1 to -14, and £ = 2 and 3 low-order (« = -1 to -5), and high-order {n = -16 
to -20) modes. We then followed the variation in frequency of each mode of degree 
£, azimuthal order m, and radial order n by slowly increasing the rotation rate, step by 
step, from to 0.1 Q.k- For a de t ailed description of the mode-following technique, we 
refer the readers to Ballot et al. (I2010h . 



3. Validity domain of perturbative methods 

The effects of rotation on the oscillation modes can be treated as a perturbation where 

the rotation rate is the small parameter. A Ist-order co rrection has been proposed by 

lL"edoux U951.), 2nd-order by ,Saio (1981). lDziembowski & Goode (1992). and Suarez et al 



(l2006h . and 3rd-order terms have been developed by ISoufi et al.1 (119981) . Within a per- 
turbative approach, frequencies are developed as 

= + ""^"'^^ + "/ + '"'^ + '"(^«/ + '"'^n/)^^ + ^(^') (9) 

The bar denotes the normalization co = co/Q^ and O = O/O^ where = (GM/Rp)^^^ 
(Rp is the polar radius). We use this normalisation since the polar radius is expected 
to be a slowly varying function of O in real stars, as opposed to /?eq. We compute the 
perturbative coefficient s C^/, 5 ^ ^, and Tl ^ from our 2D computations as described in 
lReese etal.1 (l2006h and IB allot etaP (l201o() . 



From these coefficients we calculate mode frequencies with the 1st to 3rd-order 
perturbative approximations for rotation rates ranging from O = to O.VQk and com- 
pared them to complete computations. To define the domains of validity of perturbative 
approaches, we fix the maximal departure, dto, allowed between the perturbed frequen- 
cies, oj'^'^i.ljj, and the exact ones, (lJn,{,m- For each mode and each approximation order, we 

define the domain of validity [0, Q.^], such that VQ. < Q.^ \of^^^^^JyQ) - cl)n/^m{^)\ < So). 
The precision of the observed frequencies, 6v, is related to the normalized error, 6a), 
through 6(1) = 6v/vg with = Qk/(2^)- We thus display in Fig. [T] the domains of va- 
lidity of the perturbative approximations for a frequency precision 6v = OA yuHz (spec- 
tral resolution after a hundred days) and for two types of stars with different dynamical 
frequencies: a typical y Dor star (v^ = 61 //Hz), and a typical B star (vg = 11 juHz). No- 
tice that these two plots can also be seen as the domains of validity of the perturbative 
approximations for a given star (e.g. , a y Dor) star at two frequency precisions (e.g. , 0. 1 
and 0.5 juHz). The domains of validity obviously extend to higher rotation rates in the 
bottom panel. 

We observe distinct behaviours in the high- and low-frequency ranges. In the 
high-frequency range, 2nd-order perturbative methods give satisfactory results up to 
~100 kms"^ for y Dor stars and up to ~150 kms~' for B stars. The 3rd-order terms 
improve the results and increase the domains of validity by a few tens of km s"' . These 
results are to be contrasted with those found for p modes where the domains of va- 
lidity are restricted to lower rotation rates (for 6 Set stars, which are similar to y Dor, 
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Figure 1 . Evolution of the frequencies of ^ - 1,2,3 modes, computed in the co- 
rotating frame. Perturbative approximations have been tested for a typical y Dor ( top 
panel) and for a B star (bottom panel). The green/red/blue parts of curves indicate 
that lst/2nd/3rd order is sufficient to reproduce complete calculations within an error 
6v - Q.l yuHz. Error bars on the right-hand side of each panel show 5v and 10 x 5v. 
The magenta lines indicate lo - 2Q.. For each plot, the bottom x-axis and left y-axis 
show dimensionless units, whereas the top x-axis and right y-axis show physical 
units. 
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Reese et al. (120061) find ~50-70 km s as a limit for perturbative methods) and where 



3rd-order terms bring little improvement. The rather good performance of perturba- 
tive methods at describing high-frequency g modes indicates in particular that the 2nd- 
order term gives a reasonable description of the centrifugal distortion. This might be 
surprising considering the significant distortion of the stellar surface (/?eq/Pp - 1-08 at 
O = 0.4Qk)- Actually, the energy of g modes is concentrated in the inner part of the 
star where the deviations from sphericity remain small. As a result, g modes "detect" a 
much weaker distortion than p modes, thereby making them amenable to a perturbative 
description. A particular feature that induces a strong deviation from the perturba- 
tive method concerns mixed pressure-gravity modes that arise as a consequence of the 
centrifugal modification of the stellar structure. For example, we found that, above a 
certain rotation rate, the £ = 3,n = -I mode becomes a mixed mode with a p-mode 
character in the outer low-latitude region associated with a drop in the Brunt- Vaisala 
frequency. 

However at low frequency, the domains of validity of perturbative methods are 
strongly reduced: for y Dor stars, 2nd-order perturbative methods are only valid below 
~50 kms^^. Indeed, perturbative methods fail to recover the correct frequencies in 
the inertial regime co < 2Q. (delimited by a magenta curve). In particular, we observe 
that increasing the tolerance da> between the top and bottom panels brings very little 
improvement in the inertial domain. As detailed in Sect. 15. 1[ this can be attributed to a 
change in the nature of modes in this regime. 



4. Regularities: Period spacing 

For a nom-otating star, the period spacing AP,,/ = Pn+i,e,o - P,i/,o> where Pn/.m - 
27r/w„/^m> is known to be asymptotically independent of n, in the absence of strong 

structure gradients. More precisely, AP„ / ~ APf = Pp/ s /fji + I), where Po = 2n^ (J Ng/r 
only depends on the structure of the star ( Tassoull[l980h . For a rotating star this is not 



true. AP„/ „j also depends on Q, n, and m. Nevertheless, for a given m, by considering 
sufficiently high values of n{n> 10) we realise that the period spacing depends mainly 
on a reduced parameter, rj = 2Q./a>. Results for f = 1 are displayed in Fig.|31 

This parameter ij appears naturally within the so-called traditional approximation 
(TA). Commonly used in geophysics dEckart I960), it consists in assuming a spherical 
symmetry for the star and neglecting the tangential component of the rotation vector 
in Eq. |6] By doing this, the problem becomes separable in the radial and latitudi- 
nal coordinates. The equation governing the radial dependence of the mode is similar 
to the nonrotating problem with {{£ + 1) terms replaced by /lf,„,(?7), where /t^m are 
the eigenvalues o f the Laplace's tidal eigenvalue problem in latitudinal direction {e.g., 
lUnno e t al.lll989h. The TA has been used for computing g-mode frequencies in stars 
(e.g.,!Berthomieu et al. 1978; Lee & Saio 1987, 1997; T ownsend 2003). Within the TA, 
for sufficiently high values of n, Pn,e,m ~ nPol ^J^^,mi1l)■ We then derive the spacings 



Po 



AP ^ AP(77) « — . (10) 



The star is assumed to be spherical for the TA but spheroidal for the complete computa- 
tions. Thus, in order to compare the results, we decided to consider that the polar radius 
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Figure 2. Period spacings for { - I modes as a function of 77. Pluses correspond 
to complete computations, dots to computations on a spherical model with a full 
treatment of the Coriolis force. Modes with n = -10 to -14 are used. Solid lines 
correspond to the asymptotic relation obtained within the TA. 



of the distorted models is equal to the radius of the spherical model. The relations (fTOl) 
are plotted in Fig. |2] for £ = I. They agree very well with our complete computations. 
We only see a departure between the asymptotic TA and exact frequencies for m - -€ 
modes for high values of 77. For ( = 2 and 3, we get very similar results (not shown). To 
investigate the origin of this discrepancy we performed complementary computations. 
We computed ( = \ g-mode frequencies with TOP, without any simplification to the 
Coriolis term, but using a spherical model. The resulting period spacings are plotted 
with dots in Fig.|2] They agree well with the TA even for m = which demonstrates 
that the remaining difference originates in the centrifugal distortion of the 2D models. 



5. Spatial structure of the eigenmodes 
5.1. Inertial domain o) < 2Q. (77 > 1) 

In Sect. [3l we suggest that the failure of the perturbative method in the inertial regime 
is due to a change in the mode nature. In this regime, we observed changes in the mode 
cavity that are apparently not taken into account by the perturbative method. Modes in 
the inertial regime do not explore the polar region and the angular size of this forbidden 
region increases with 77. This is illustrated in Fig.[3]for a particular mode. Such a drastic 
change in the shape of the resonant cavity has a direct impact on the associated mode 
frequency. As perturbative methods ignore this effect, they cannot provide an accurate 
approximation of the frequencies in this regime. 

This interpreta tion is supported by t he an alytical expression of the forbidden re- 



gion determined by lDintrans & RieutordI (12000 ) for gravito-inertial modes. Indeed, for 



frequencies oj < 2Q, the Coriolis force becomes a restoring force, and modes become 
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Figure 3. (Left) Meridional distribution of kinetic energy of the g mode (f - 
3,m = -l,n - -18) in a nonrotating star (Ri^ht) The same forQ = 0.7Qk- Magenta 
lines indicate the critical sur face r = 0. See also lBallot etall (l2010h . 



mixed gravity-inertial modes. With a spherical model, and within the anelastic and 
Cowling approximations, they have shown that gravito-inertial waves only propagate 
in the region where 

r = r'^oj^lNl + (inf - 0?] - (IQNozf > 0. (1 1) 

This implies that, when 77 > 1, a critical latitude 6c = arcsin(l/?7) appears above which 
waves cannot propagate. Even though this expression does not strictly apply to our 
nonspherical geometry, we have superimposed the critical surfaces, F = 0, over the 
energy distributions of our eigenmodes (Fig.O. Without rotation, there is only a small 
circle close to the centre, corresponding to the classical turning point oj - No- For the 
mode with a» < 2Q, the polar forbidden region delineated by F = agrees pretty well 
with the energy distribution of the mode. 

In real stars, which have convective cores, the gravito-inertial nature of modes 
could become even more crucial since the gravito-inertial waves propagate in convec- 
tive regions when g modes cannot. 

5.2. Rosette modes 

Even if they are not described by a unique spherical harmonics, modes with 7/ < 1 
are generally strongly dominated by one degree. Hence, they look like "distorted" 
spherical harmonics (of course, avoided crossings sometimes create a strong coupling 
between modes leading to more complicated patterns). However, for a few modes, 
such as (f = 3,m = 0,n = -3), the mode structure changes quickly as the rotation 
rate increases (for this specific mode, changes are important at rotation rates as low as 
O.IQk): the energy distribution focuses around a "rosette" pattern. Figure |4] displays 
such a rosette mode. 

Following the same approach as lLignieres & GeorgeotI ( 20081 12009) for p modes. 



Prat et al. (in preparation) have developed a ray theory for gravito-inertial modes in 
rotating stars. Using these developments, we computed the ray dynamics for the fre- 
quency of the {£ = 3,m = 0,n = -3) mode at Q = 0.38nK ((^ - 1-V2). The Poincare 
surface of section (PSS) is plotted in Fig. |4] (for more details on the PSS and notions 
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Figure 4. (Left) Meridional distribution of kinetic energy of the g mode ({ = 
3,m - Q,n - -3) for Q. = O.SSOk- The red Une is a stable periodic orbit obtained 
from ray theory. (Right) PSS corresponding to the same frequency (w - 1.72). The 
PSS is computed at - njl (equatorial plane), kr denotes the radial component of 
the wave vector. 



related to ray dynamics, see iLignieres & Georgeoj 120091 ') . We first see an unexplored 
region in the core (r < 0.08), corresponding to the region where o) < No- The most 
striking features of this PSS are the two large islands formed around stable periodic 
orbits. One of these two stable orbits is superimposed over the mode distribution in 
Fig- HI and we see that the energy of the mode is distributed around it. By exploring the 
spectrum with TOP in the vicinity of this mode (i.e. m = modes with a> ~ 1.72), we 
discover other rosette modes the energy of which is also distributed around the same 
rays, but with more nodes. 

Thanks to ray theory, we also discover another family of rosette modes, around 
CO = 1.18, corresponding to 2-period islands, that are also found with TOP. All of these 
modes still exist even when the centrifugal distortion is neglected, i.e. within a spherical 
approximation. 
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